Gdl 2 : A New Ferromagnetic Excitonic Solid? 
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The two-dimensional, colossal magnetoresistive system GdI-2, develops an unusual metallic state 
below its ferromagnetic transition and becomes insulating at low temperatures. We argue that 
this geometrically frustrated, correlated poor metal is a possible candidate for a ferromagnetic 
excitonic liquid. The renormalized Fermi surface supports a further breaking of symmetry to a 
charge ordered, excitonic solid ground state at lower temperatures via order by disorder mechanism. 
Several experimental predictions are made to investigate this unique orbitally correlated ground 
state. 

PACS numbers: 71.45.Lr, 71.30.+h, 75.50.Cc 



c3 



i 

C 
O 

o 



> 

O 

CO 

o 

oo 
o 



% 



Correlated electronic systems with geometric frustra- 
tion have spawned intense interest recently, because of 
the startling diversity of physical properties they man- 
ifest. The crucial roles of electronic correlation and 
frustration, working in tandem, are nevertheless poorly 
understood from a fundamental viewpoint: strong cor- 
relations generically drive Mott transitions (accompa- 
nied by orbital/magnetic order) and manifest spectac- 
ular changes in physical responses under minute pertur- 
bations [1]. In correlated systems, geometric frustration 
gives rise to large, exponential degeneracy of classical or- 
dered states, and to unconventional order via an order 
by disorder mechanism [2| . The variety of new, hitherto 
unexplored, phases of matter emerging in systems where 
both correlation and frustration are operative is a fasci- 
nating, more-or-less open problem. 

The layered (in XY plane) ferromagnetic (FM) system 
Gdh, studied by several authors [3, U, la, |°|, may turn 
out to be an interesting example. Gd ions in the 4/ 7 5d 1 
state form a triangular lattice, separated by iodine lay- 
ers. Itinerant electronic states in the partially filled ri- 
bands coupled to localized /-electrons drive an FM state 
(T c ~ 300 K). Interestingly, Gdh shows significant mag- 
netoresistance [a] (40% per Tesla) close to T c [lflj. In 
addition, it is iso-structural to the well known dichalco- 



genides with hexagonal layered structure 11|, showing 
charge density wave (CDW) order at low temperatures. 
In fact, Gdli is a bad metal below T c and becomes in- 
sulating for T < 80 K. Might this be an indication that 
the features seen in dichalcogenides also occur in Gdl%l 
The unusually high (significantly higher than p^f^t) 
resistivity of Gdli even in the metallic state is strongly 
at odds with a simple FM s(d) — f lattice model, where 
FM order implies a good metal for a translationally in- 
variant system. In fact, experiments yield kpl < 1 (I 
is electron mean- free path) in Gdli, indicating a strong 
inelastic scattering of carriers and resulting in an incoher- 
ent metallic state, without Fermi liquid quasi-particles. 



On the other hand, the FM spin fluctuations |l2j are well- 
described by a classical picture of two-dimensional (2D) 
FM Heisenberg model. The saturation magnetic moment 
in the FM state is about 7.33fiB per Gd atom, consider- 
ably lower than the maximum attainable value (8/is), in- 
dicating the relevance of the competition between strong 
correlation effects and itinerancy among the 5d electrons 
of Gd. The above observations force one to seek addi- 
tional, non-magnetic, strong correlation effects deep in 
the FM state to reconcile magnetism with transport in 
Gdl2- Moreover, owing to its 2D structure, the role of 
correlation is expected to be very important in Gdl2- 
Strong magnetic fluctuations found experimentally |l3| 
and a series of magnetic ground states obtained in a cor- 
related electronic model |P|, also underline the relevance 
of Gd d-shell correlation. 

Band structure results for Gdl2 |3, |9( indicate a spin 
splitting of the conduction band. The crystal field in 
Gdl2 is relatively small and the 5c? orbitals are weakly 
split and partially occupied. Emergence of an insulat- 
ing ground state in with three partially filled d-bands 
is surprising, since the band-filling is off-commensurate 
(less than 1/3 filling). This suggests that the T — > 
insulator must be associated with an additional broken 
symmetry. The high-T bad-metal phase implies that this 
must emerge from an instability of an incoherent metallic 
state, pointing to the strong correlation limit. This raises 
deeper questions: (i) what is the specific nature of the in- 
clastic scattering processes leading to the "bad metal"? 
(ii) How does such an incoherent FM state evolve into a 
T — * insulator? Is any additional symmetry breaking 
(in charge, orbital sectors) involved? 

These observations should go hand-in-hand with the 
evolution of the correlated electronic structure of Gdl2 
with temperature. In this letter, we address these issues 
in detail. Using dynamical mean-field theory (DMFT), 
we solve a multi-orbital Hubbard model with the LSDA 
band-structure of Gdl2 [9]. We show that strong cor- 
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FIG. 1: (Color online) The LSDA+U DOS of the five 
Gdd— bands in the majority and minority spin sector for 
(7 = 6eV 



relations localize two degenerate bands with d x 2_ y i 1 d xy 
orbital character, leaving a renormalized band with d z 2 
character to cross Ep . The renormalized Fermi surface is 
strongly modified by correlations, eroding the LSDA hole 
pockets (see below). Further, we argue that the shape of 
the correlated Fermi surface naturally drives a charge- 
ordered (CO) ground state with \/3 x \/i order on the 
triangular lattice (breaking the underlying translational 
symmetry). Using LSDA+DMFT, we propose a simple 
effective model that allows us to suggest experimental 
signatures of such a ground state. 

The LSDA band structures obtained using LMTO 
method [14] shows that the most relevant bands crossing 
Ep are mainly derived from the bd orbitals (occupations 
of d x 2_ y 2, d xy , d z 2 are 0.31, 0.31, 0.30 respectively), split 
about 0.9 eV into majority and minority spin sub-bands 
by strong exchange interaction (Jfd) with the completely 
spin-polarized 4/ shell. The corresponding hexagonal 
Fermi surface, exhibits six hole pockets (not shown). 

LSDA+U calculations were performed treating Gd / 
states as completely spin-polarized quasi-core states (U 
applied only to Gd d states). For the double counting 
term the so-called "atomic limit" was used [6j, [Z| . We find 
(Fig. [T]) that d x 2_ y 2, d xy and d z 2 bands are pulled down 
by static, orbital correlations, and the d xy , d x 2_ y 2 bands 
are localized. However, LSDA+U overestimates band- 
splittings, while LSDA cannot split partially occupied 
bands. A reliable picture requires a full incorporation of 
dynamical d-band correlations, using LSDA+DMFT, to 
which we turn below. 

The correlated multi-orbital model Hamiltonian for 
Gdli is described by H = H + H int , where 
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describes the correlation part. In addition to U, U' act- 
ing in the d-manifold, we include Jfd for the scatter- 
ing and spin polarization of the d-bands by well-localized 
Gd 4 /-states (treated classically, i.e., neglecting spin- flip, 
as S = 7/2). As parameters relevant to Gdl2, we choose 
U = 7.0 eV, U' = 5.0 eV, and J fd = 1.5 eV, with the 
LSDA bandwidths between 5 — 6 eV. In addition to scat- 
tering induced by U, U', Jfd acts like a classical scat- 
tering potential on the d-band states. We consider the 
FM phase without additional (charge, orbital) symme- 
try breaking and use multi-orbital iterated perturbation 
theory (MO-IPT) [15j to solve the multi-orbital impurity 
problem in the d-sector. LSDA+DMFT has been used 
with good quantitative success in many 3d-oxides |15l.ll6j|. 
The strong scattering induced by Jfd is not accounted for 
in multi-orbital DMFT and is therefore combined with 
the solution of the Falicov-Kimball (FK) model in the 
local approximation 16]. 

DMFT rcnormalizes the LSDA results in two steps. 
First, the multi-orbital Hartree self-energy renormal- 
izes the relative (LSDA) band positions depending upon 
their occupations and (energy) separations. More im- 
portantly, the frequency-dependent self-energy causes 
spectral-weight transfer across large energy scales, drasti- 
cally modifying LSDA spectra. Its importance is directly 
seen in photoemission experiments on a host of correlated 
systems (l7|: DMFT generally gives good quantitative 
agreement with photoemission. While LDA(LDA+U) 
generically gives correct ground states with orbital and 
magnetic order for weakly correlated(correlated) solids, 
its inability to describe both the narrowed quasipar- 
ticle bands and high-energy satellite features is well- 
documented 17 1 . 
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Our DMFT results are drastically different: 
Hartree self-energy pushes the degenerate d x 2_ y 2, uxy 
bands down, and for U' = 6 eV, moves them completely 
below Ep, leaving only the narrowed d z 2 band to cross 
F,p (Fig. [2|) . Smaller values of U' did not yield this selec- 
tive localization. As alluded to above, large scale spectral 
weight transfer to the Hubbard bands drastically modi- 
fies the LSDA DOS. In contrast to multi-orbital systems 
showing correlated Fermi-liquid behavior at low-T, no 
sharp quasiparticle peak (i.e., /mE(w — ► 0, T — > 0) + 1 0) 
is obtained in the DOS here. This implies [18( an inco- 
herent metallic behavior, with lifetimes so short that the 
quasiparticle concept loses meaning. This agrees with the 
high metallic resistivity in Gdl^, which can be classified 
now as an orbital selective "bad metal" . 

In fact, with Mott localization of the d x 2_ y 2, d xy 
bands, the problem effectively reduces to a FK model 
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FIG. 2: (Color online) The spectral density for d xy and 
d x 2_ y 2 derived bands from DMFT calculations (lower panel). 
The upper panel represents the same for d z 2 derived bands. 



at low energy where the itinerant d z i electrons scat- 
ter off the localized d x 2_ y 2 and d xy electrons. The in- 
frared X-ray edge singularities rigorously known to ex- 
ist for the FKM in D = 



181 ] give an infrared diver- 



gent local, inter-orbital excitonic susceptibility, \ab (w) — 
/^<r[aL(t)^(t);fol(0)a iCT (0))e- t , implying the exis- 
tence of soft, inter-orbital (a, (&) = d z 2, (d x 2_ y 2, d xy )) ex- 
citonic modes. Coupling itinerant d z 2 electrons to these 
soft modes then destroys Fermi-liquid coherence leading 
to a "bad metal" . Given finite T = entropy in this 
bad metal, we expect that a broken symmetry phase will 
preempt this unstable state at lower T. The quasi-nested 
regions of the renormalised Fermi surface (Fig. [3]) indi- 
cate that this broken symmetry state is likely to be a CO 
state. 

We can now draw the following straightforward con- 
clusions. Since the d x 2_ v 2, d xy bands are pushed below 
Ep by multi-orbital DMFT, the small hole pockets found 
in the LSDA disappear in the renormalized Fermi sur- 
face, which now corresponds to the correlation-narrowed 
d z 2 band crossing Ep . The DMFT orbital occupations 
are altered to 0.35,0.35,0.22 for the d x 2_ y 2 1 d^and d z 2 
bands. Similar features are seen in DMFT work [19j on 
Na x Co02, in agreement with ARPES [20j showing no 
hole pockets. We predict that (i) an ARPES measure- 
ment performed for 80 K < T < T c would show up a 
single hexagonal Fermi surface sheet, and (ii) ARPES 
lincshapes will be anomalously broad without any Fermi- 
liquid peaks at low energy, reflecting the "bad metal" 
state. 

A very interesting aspect of the renormalized Fermi 
surface, however, is that, similar to dichalcogenides 11|. 




FIG. 3: (Color online) The dispersion of the LSDA+DMFT 
band crossing the fermi level and the corresponding fermi sur- 
face are shown. 



it exhibits a built-in tendency to favor emergence of low- 
T CO state of the v^ x 



111 ]. The movement 



\/3 type 

of the van-Hovc-like peak in the d z 2 DOS very close to 
E F in the LSDA+DMFT (in stark contrast to LSDA 
and LSDA+U results) further corroborates this argu- 
ment. Such features in Fermi surface and DOS may lead 
to unconventional CO states with gapless excitations [8(] . 
Rigorous results for the existence of such CO states exist 
in the context of FK model (see below). 

The essential results can be understood in a simpler, 
effective model for the FM state at low energy. With elec- 
trons in d x 2_ y 2, d xy bands unable to hop, we introduce 
an effective Hamiltonian, defined by 



H eff = ^2 € z 2 (^)dl2j,d z 2 tk + A^(n M 2 -^n ijQ ) 

k i a 

+ U\ J^ hi, a h^ z 2 (3) 



where a=d x 2_ y 2, d xy . We solve H e ff without fur- 
ther approximation within DMFT, with the LSDA oc- 
cupations (0.31, 0.31, 0.30 for d xy , d x 2_ y 2 and d z 2 or- 
bitals). The inter-orbital correlation U\ reduces the DOS 
at Ep considerably (Fig. 2]) and a bad metallic phase, 
with Ep in the low-energy pseudogap emerged. The 
DOS shows a gap for U\ = 0.70 eV; the Fermi level, 
though, is never in the gap owing to the partial occu- 
pancy of the band. The local inter-orbital susceptibilities, 
Xz 2 ,xy "(u),x z 2 tX 2_ y 2(u) exhibit power-law divergences 
(in the D = oo FKM [l8j), yielding incoherent metallic 
behavior. 

Reduction to an effective FKM has further, concrete 
consequences. Since only the d z 2 states are itinerant, the 
local constraint of single occupancy on each site in -ff e // 
leads to the emergence of varied, band-filling dependent 
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FIG. 4: (Color online) The density of states for the effective 
band from the model Hamiltonian Eqn. (3). The red and 
green curves are for U\ — 0.70 and 0.30eV. The blue one is 
the DOS for the tight binding fit to the LSDA band that 
crosses the Fermi level. The chemical potentials are set at 



CO states. This is rigorously known [2l| for FKM, and, 
with three (orbital) components in Hamiltonian Eqn. (3), 
we expect related CO patterns to unfold as a function of 
band filling. Indications of this may already have been 
seen in Gdli [22j. We propose that such an "ordered" 
state also drives Gdl 2 insulating at low T. Moreover, the 
detailed shape of our LSDA+DMFT Fermi surface sug- 
gests an in-built propensity towards an unconventional 
CO state with an anisotropic gap at Ep. 

The Mott-Hubbard localization of d x 2_ y 2, d xy band 
states found above leads to a local U(\) invariance forc- 
ing [13 (d\ 2 d a 



identically (a=d x 



-y^, u<xy 



d xy ). This 



implies rigorous absence of inter-orbital excitonic order 
in the bad metal. Technically, this also precludes use of 
usual Hartree-Fock RPA approach for broken symmetry 
phases here. Interestingly enough, exactly this circum- 
stance generates [18[ infrared singular, local, excitonic 
modes in the uniform phase, as shown above indicating 
that the FM bad metal in Gdl 2 is an excitonic liquid. 
The LSDA+DMFT results do show, however, that the 
uniform phase is intrinsically unstable (given soft exci- 
tons [24[) to an inter-orbital excitonic order at low tem- 
perature. We propose that the unconventional CO state 
should also be viewed as an unconventional excitonic 
solid. We speculate that this could turn out to be a new, 
specific microscopic manifestation of a seemingly more 
generic feature in dichalcogenides ll|. Suitable external 
tuning parameters (e.g., pressure, doping) may melt this 
state into an excitonic liquid. With discrete (Ising) na- 
ture of orbitals, frustration on a triangular lattice may 
lead to partially ordered solid-like, or more exotic dimer 
or plaquette ordered phases before the unconventional 



state melts completely. 

The structural similarity of Gdl2 to other dichalco- 
genides like 2H-TaS2, 2H-TaSe2, etc, suggests related 
phenomena in Gdl 2 . In particular, it is possible that 
the unconventional CDW or excitonic solid (or liquid) 
states manifest themselves - favorable cases would be 
those that are near the metal-insulator transition [25| . 
High pressure or chemical doping could then tune the sys- 
tem into a very anomalous metallic state with strong, un- 
conventional, excitonic liquid correlations. Whether such 
excitonic-liquid correlations also drive unconventional su- 
perconductivity (USC) [ll| | is an additional, fascinating 
issue in this context. 

To conclude, based on LSDA+DMFT calculations, we 
propose that the FM "metallic" phase of Gdl 2 is an ex- 
citonic liquid. At low T, this is argued to order into an 
unconventional excitonic CDW insulator. Our approach 
can also be used fruitfully for the investigation of other 
dichalcogenides of great current interest. 

LC thanks the Emmy-Noether Program of the German 
DFG for support. MSL thanks the MPIPKS, Dresden for 
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